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HIGHLIGHTS 


• The effect of particle shapes on pulverized fuel particles was analyzed. 

• Theoretical considerations were tested and applied on experimental data. 

• The particle shape was found to have an effect on the burning rate of fuel particles. 

• Current models should be adapted for non-spherical particles. 
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Biomass combustion becomes increasingly important in nowadays energy production. One of the major 
differences between biomass and coal is the particle shape. Therefore the effect of particle shape varia¬ 
tions was calculated to determine the effect of spheroidal or cylindrical particle shapes when determin¬ 
ing char burning kinetics of pulverized fuel particles from pyrometric measurements. A model frequently 
used for char kinetics determination was adapted to spheroidal and cylindrical particles under typical 
combustion conditions. The impact of the particle shape on the energy balance, which is needed for 
the determination of kinetics parameters, was calculated for typical particle sizes and length to diameter 
ratios. As the calculations show, the influence of particle shape effects is not negligible for the determi¬ 
nation of biomass char burning parameters by pyrometric measurements. The application of the adapted 
model on measured data shows, that particle shape effects can be resolved by pyrometric techniques. The 
results presented indicate that a sound analysis of the particle burning rate should include a three- 
dimensional particle shape measurement. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Due to the demand for carbon-neutral fuels, the amount of 
biomass replacing coal in pulverized fuel boilers is increasing [1]. 
In pulverized fuel (pf) combustion, the combustion of char (fixed 
carbon) is the time-determining step. Design of pf-boilers, which 
is usually supported by CFD-simulations, requires an exact predic¬ 
tion of the char burning rate, both for the determination of spatial 
heat release and for sufficient burn out. 

At least three competing techniques to determine char burning 
kinetics are known in literature. TGA experiments have been 
performed on a variety of fuels (coals, raw and torrefied biomass) 
under air-fired and oxy-fuel conditions 2-4]. In TGA experiments, 
the temperature is usually kept low enough to assure chemically 
controlled (Zone-I) combustion conditions, but at low heating rates 
which are not typical for pf-combustion. Typical heating rates are 
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easily achieved in drop tube reactors (DTR) or entrained flow reac¬ 
tors (EFR), which have been applied to study the combustion and 
gasification rates of coal and biomass char under air-fired and 
oxy-fuel conditions [2,5 . In DTRs the fuel particles follow a linear 
trajectory at a nearly constant velocity, which simplifies the deter¬ 
mination of residence times. A common technique to determine 
the reaction kinetics parameters is sampling of solid combustion 
residues. These samples are then used to measure the degree of 
conversion by the so-called ash tracer method [6]. 

1.1. Optical determination of burning rates: state of the art 

Another method to determine the apparent reaction kinetics of 
char particles in DTRs is two-color pyrometry. This method is 
based on measuring temperature and diameter of a single burning 
char particle under well-defined conditions. From these experi¬ 
mental data, the energy balance around the particle is calculated. 

Several research groups have deployed ratio pyrometers in two 
or three-color versions to determine char burning rates. At Sandia 
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Nomenclature 




pre-exponential factor (mol/(m 2 s)) 

Qconv 

convective heat flux R (W/m 2 ) 

A 

area (m 2 ) 

Q inert 

heat flux into the non-reacting particle (W/m 2 ) 

A C yl 

cylinder surface (m 2 ) 

Qi 

heat flux (W/m 2 ) 

A s 

outer surface area (m 2 ) 

Qrad 

radiative heat flux (W/m 2 ) 

c p 

heat capacity of char (J/(kg K)) 

Qreact 

surface-specific heat of reaction (W/m 2 ) 

d 

minor particles’ axis (ellipsoid), diameter (cylinder) 

Qi 

energy of type i (J) 

d A 

surface diameter (pm) 

R 

molar gas constant (J/(mol K)) 

dchar 

characteristic diameter (pm) 

Spx 

pixel size (pm) 

d P 

particle diameter (pm) 

t 

time (s) 

dsphere 

diameter of particle with equal projection area (pm) 

Tf 

film temperature (K) 

d v 

volume diameter (pm) 

T g 

gas temperature (K) 

E a 

activation energy (J/mol) 

T P 

particle temperature (K) 

h 

heat transfer coefficient (W/(m 2 K)) 

T w 

wall temperature (K) 

AH 

heat of reaction (J/mol) 

v p 

particle velocity (m/s) 

L 

largest dimension of particle (m) 

V C yl 

cylinder volume (m 3 ) 

n 

reaction order (-) 

v P 

particle volume (m 3 ) 

tlpx 

number of pixels 

X 

variable (-) 

Nu 

Nusselt number (-) 

z 

distance to fuel inlet (EFR) (m) 

Po 2 ,s 

0 2 -pressure at a particles’ surface (Pa) 

a 

heat transfer coefficient (W/(m 2 K)) 

Q 

burning rate (mol/(m 2 s)) 

£ 

emissivity (-) 

Qcyl 

burning rate assuming cylindrical particle (mol/(m 2 s)) 

2 

thermal conductivity (W/(m K)) 

Qell 

burning rate assuming elliptical particle (mol/(m 2 s)) 

Pp 

char density (kg/m 3 ) 

Qsphere 

burning rate assuming spherical particles (mol/(m 2 s)) 

o 

Stefan-Boltzmann constant (5.67 * 10“ 8 (W/(m 2 K 4 ))) 


National Laboratories, two-color pyrometry of burning char parti¬ 
cles in an EFR has been carried out for a long period. From the late 
1970s on, coal and biomass chars have been investigated for their 
burning behaviour [7-13]. This system is based on a coded aper¬ 
ture for particle size detection. As the particle size is deduced from 
intensities measured by photo multipliers (PMT), the particle 
shape is assumed to be spherical. Another system based on size 
measurement by a coded aperture has been developed by Cope 
et al. [14], also using PMTs to measure the intensity of burning coal 
char particles. Other groups developed systems using fibre optics 
to guide particle radiation to PMTs, e.g. Spliethoff et al. 15] and 
Joutsenoja et al. [16,17], both using pyrometric data from single 
burning particles to determine char burning rates. All pyrometers 
mentioned so far have in common to measure in radial direction 
perpendicular to the coal streak in an EFR or DTR. Therefore only 
a momentary of each particle passing the pyrometers field of view 
is detected. Typically, a number of particles sufficient for reliable 
statistics are collected at a specified distance to the solid fuel inlet, 
and several measurements at different distances (residence times) 
are required to cover the burning history of a fuel sample. Another 
approach has been chosen in the work of Levendis and his group. 
The pyrometer system, described in detail in [18], measures along 
the axis of the EFR, which allows to measure the complete trajec¬ 
tory of a single particle. Thus, beneath the determination of particle 
temperatures [19], the burning time can be measured directly 
[20,21]. Particle burning rates have been derived from these mea¬ 
surements [20], and literature data has been tested to fit the mea¬ 
sured burning times and particle temperatures of different coal 
samples with good agreement [21]. 

A system based on intensified (I)CCD cameras has been used to 
measure particle temperature and size by two-color pyrometry by 
Hackert et al. [22]. As the ICCD imaging technique provides spa¬ 
tially resolved 2d-images of burning particles with a resolution 
of 10 pm or less, particles of typical size for pf-firing (d p = 50- 
100 pm) can be resolved spatially. The pyrometer system has been 
used for further investigation of coal and biomass chars, an exam¬ 
ple on particle temperature measurement of coal chars in oxyfuel 
atmospheres is given in [23]. 


While coal particles in pf-applications are typically compact and 
nearly spherical, even for those particles shape effects have been 
considered in literature [24]. Biomass particles originating from 
woody and herbaceous fuels show severe deviations from the spher¬ 
ical shape [25-27]. As will be described in the following, the particle 
shape has to be considered when deriving char burning kinetics 
from pyrometric measurements, as the energy balance, which is 
the key element of burning rate calculations, is affected by the par¬ 
ticle shape. As the fibrous structure of biomass particles leads to par¬ 
ticle shapes which can be assumed to be cylindrical or ellipsoidal, 
both shapes are tested in theory and on exemplary measurements. 


2. Char burning kinetics from pyrometry 

2 A. Energy balance around burning particles 


The key element of burning rate determination by pyrometry is 
the energy flow balance at the outer surface of a particle: 


Q react — Qrad T - Qconv 


"T Qtnert 


m 2 s 


0 ) 


As shown in Eq. (1), the energy produced by reaction is 
balanced by the radiative heat transfer between the particle and 
its surrounding, the convective heat transfer between particle 
and gas phase and the temperature change of the inert particle 
(the non-reacting char and ash). In the following, the participating 
energy fluxes are calculated for particles of spherical, prolate 
ellipsoidal and cylindrical shape (Fig. 1). 

In general it is assumed, that the surface temperature is 
constant and the heat and mass transfer are homogeneous, so that 
the surface dependency of q, = dQ.//(dt * dA) can be expressed by 
setting dA = A s , where A s denotes the particles outer surface. Under 
the given assumptions, the expressions for the reaction energy and 
the radiative heat transfer are shape-independent, as shown in 
Table 1. The reaction part of the heat flux is calculated from the 
surface reaction rate q and the heat of reaction AH: 

Qreact — QAH. (2) 
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Fig. 1 . Shapes of the model particles considered, L/d = 2.5. 


Table 1 

Volume/surface ratio for spherical, prolate ellipsoidal and cylindrical particles. 



Sphere 

Prolate ellipsoid 

Cylinder 

V P IAs (m) 

d 

6 

Ld 

Ld 

2(d+2L) 


The radiative heat transfer is calculated from the Stefan-Boltz- 
mann constant cr, the particles emissivity e, its temperature und 
the wall temperature (possible effects of gas species affecting the 
radiative heat flux are usually neglected): 

q r ad = <7s(T 4 p -T 4 w ). (3) 

To calculate the convective heat transfer ( q con v = h(T p - T g )) 
between the particle and the gas phase, the heat transfer coeffi¬ 
cient h can be calculated as: 


where Nu is the Nusselt number, X is the average thermal conduc¬ 
tivity of the boundary layer, and d C har is the characteristic diameter 
of the particle. For pf-particles in EFR reactors the assumption of 
Nu = 2 is valid as the relative velocity between particle and gas 
phase is usually negligible [9]. 

The calculation of the energy balance presented above is based 
on the single-film model, meaning that direct oxidation of the par¬ 
ticle surface occurs, where CO and/or C0 2 are produced. The CO is 
oxidized in the free stream, not affecting the particles heat balance 
[28]. This assumption is acceptable for particles in the size range of 
50-100 pm. 

The effect of Stefan flow produced by the simultaneous produc¬ 
tion of CO and C0 2 is known to influence the char burning rate by 
several percent depending on the conditions applied [29,30]. For 
particle temperatures below 2000 K it seems acceptable to neglect 
this effect leading to only slightly underdetermined burning rates, 
as the results in [29] show. 

2.2. Energy balance for non-spherical particles 

Two summands in Eq. (1) depend on the particle shape, q inert 
and q CO nv The heat flux into the non-reacting particle (remaining 
carbon and ash) is shape-dependent, as the particle is assumed 


to be of constant density resulting in the necessity to calculate 
the ratio V P IA S . To calculate q inert = ( V p /A s )p p CpV p dT p /dz , the den¬ 
sity of the char p p , its heat capacity c p , the absolute particle velocity 
v p and the temperature change depending on the distance to the 
fuel inlet dl p \dz are taken into account. The last two parameters 
are typical for optical EFR-experiments, were the absolute particle 
velocity can be measured easily and the temperature change is cal¬ 
culated from the change in the average temperature along the par¬ 
ticles trajectory. The resulting volume/surface ratios ^ are given in 
Table 1, where L and d denote the diameter and height (cylinder) 
respective long and short axes (Fig. 1). 

Shape coefficients are necessary to calculate the heat transfer 
coefficient for particles of arbitrary shape. Furthermore, d char has 
to be chosen properly to match the shape of the particle, as shown 
by Wadewitz and Specht [31 ]. The sieve diameter, which is usually 
applied to characterize the particle size of pulverized fuel, is not 
sufficient without correction. As an alternative, the surface diame¬ 
ter d A and the volume diameter d v were found as best alternatives 
when calculating the heat transfer by convection using: 

Nu = 0(Adaa L (5) 

A A 

The factor (aA)/2 is called shape coefficient. For the three parti¬ 
cle shapes considered here, the resulting shape coefficients are 
listed in Table 2. These are chosen for particles following the sur¬ 
rounding gas stream slip-free, which is a typical assumption in 
EFR-experiments. 

Based on the formulae for shape coefficients and the Nu-Corre- 
lation given in Table 2 and Eq. (5), the convective heat flux: 

Qconv — -7 (Tp — Tg) (6) 

Clchar 

is calculated to close the heat flux balance around burning particles. 
In the following, the heat balance is calculated exemplary to show 


Table 2 

Shape coefficients for quiescent particles. 


Shape 

Shape coefficient (aA)/A (m) 

Refs. 

Sphere 

2nd 

[35,36] 

Prolate ellipsoid (L > d) 

2nd-\/ (L/d) 2 -l 

[35,36] 


In (L/d+V(L/d) 2 —l) 


Cylinder (L > d) 

[8 + 6.95(L/d) a76 ]f 

[36] 
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the influence of the particle shape on burning rates when the deter¬ 
mination of burning rates is carried out by pyrometric particle tem¬ 
perature and diameter measurements. A number of assumptions 
are necessary to insert all parameters which are needed for the cal¬ 
culation. Typical values were chosen from [11,17,20,23]. q inert was 
calculated using p p = 700 kg/m 3 , c p = 2000 J/(kg K) (both justified 
in 11 ]), v p = 3 m/s and dT p \dz = 2500 K/m, both taken close to the 
upper limit in EFR-experiments [11,23] to show the maximum 
impact of these factors. From the same literature, the temperatures 
were chosen as T p = 1800 K, T g = 1600 K and T w = 500 K, the latter 
being typical for experiments in EFRs with quartz glass walls. The 
char emissivity was set to s = 0.8 [9,11]. The thermal conductivity 
of the gas phase was set to X = 0.115 W/(m K), which is also typical 
for EFR experiments. The selection of variations in particle size 
respective the ratio L/d is based on the typical properties of particle 
size measurements in pyrometry. 

2.3. Theoretical analysis of shape effects 

Three cases were investigated, starting with the assumption 
that particles with projected surfaces of 1.96*10 -9 m 2 ( d sphere = 
50 pm), 7.85 * 10“ 9 m 2 ( d sphere =100 pm) and 3.14 * 10“ 8 m 2 
( ds P here = 200 pm) were measured. The ratio L/d was varied from 1 
(1.05 for ellipsoids, as the equations for elliptical particles are 
not defined for L/d = 1) to 10 to compare the influence of the parti¬ 
cle shape. To study the influence of q rea ct = Qmd + Qcond + Qinert for 
assumed aspherical/spherical particles, which is balanced by the 
heat loss mechanisms discussed, the ratio q rea ct spherical/Qreact,sphere 
was plotted against L/d in Fig. 2. Furthermore for non-spherical 
particles it was assumed that the orientation of the particles axis 
is perpendicular to the direction of view. 

In the diagrams, the influences of the particle size and particle 
shape are clearly visible. In the case of d sphere = 50 pm, the heat flux 
balance of the ellipsoid is equal to the sphere (lower limit) and 
shows a monotonic increase. At L/d = 10, an increase of 17% relative 
to the sphere with equal projection area is found, leading to a lower 
reaction rate assuming a spherical particle. For the cylindrical par¬ 
ticle q reac t is 7% smaller at L/d = 1 and 12% larger at L/d =10, resulting 
in a smaller but still considerable influence on the reaction rate. 
With increasing particle diameter, the trend slightly changes when 
varying the L/d ratio. The ellipsoid shows little deviations from the 
sphere for L/d < 8 when d sphere =100 pm is considered, but the 
cylinder shows a clearly lower q reac t for small L/d ratios, which is 
equivalent to lower burning rates. A further increase in particle size 
to d sphere = 200 pm leads to permanently lower q rea ctS when the par¬ 
ticle shape is considered. For cylindrical particles, the effect reaches 
its maximum at L/d = 10, where q rea ct would be overestimated (note, 
that over- and underestimation is always referred to the spherical 
particle for comparison) by 14% assuming a sphere instead of a 
cylinder. For large L/d values, q rea ct of the prolate ellipsoid converges 
against the same value as the cylinder. 

The results in Fig. 2 clearly indicate the importance to consider 
the particle shape when calculating q rea ct for burning rate determi¬ 
nations. As the results show, q rea ct for small particles is underesti¬ 
mated for a cylindrical particle (L/d = 1) compared to a spherical 
one, while the elliptical particle shape shows a steady q rea ct at this 
point. It must be remembered that the calculation is based on the 
assumption of measuring the particle size/surface for a particle 
with its major axis orientated orthogonally to the direction of view. 
The particle dimensions are included in both the heating of the 
inertia as V p /A s (Table 1) and Nujd = [aA2 _1 ]/A, where the shape 
coefficient written in square brackets is taken from Table 2, which 
also depends on L and d. With L/d = x it follows: 

^cyl _ d X 

Acyl 4x + i- 1 ; 



L/d 


1,2 


cu 

*C3- 

« 1 


0,9 


G 


0.8 


d = 100 pm 


■ellipsoid/sphere 

cylinder/sphere 


4 6 

L/d 


10 



L/d 


Fig. 2. q react}aspherica i/q react , spher e for three different particle sizes, constant projection 
surface area. 


For spherical particles Eq. (7) becomes d sp here/6. When the rect¬ 
angular face of a cylindrical particle with L/d = 1 is considered, it 
follows that d = d spher e{nl4f 12 . Substituting this in Eq. (7) leads to 
Vcyi/Acyi being smaller than d sp herel 6. Additionally substituting L by 
dx leads to: 

oAl 1.59 

t a = ~t- (8) 

Compared to a spherical particle, where [a/U -1 ]/A = 2/d is used, 
this reduces the convective heat flux for a given particle when x = 1 
is applied. As is shown in Eqs. (7) and (8), for particles with small 
L/d the shape-dependent contributions differ significantly from 
those calculated for spherical particles, when cylindrical particles 
are assumed. 

The contributions of heating the inertia, q merf , and the convec¬ 
tive heat flux, q con v ,the two summands expressing the shape 
dependency, are compared in Fig. 3 (the radiative heat flux is 
constant at the particle surface for given T p and T w ). The influence 
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Fig. 3. q i ne) t and q cond for particle sizes of 50, 100 and 200 pm. 

of convective heat flux is largest for d sp here = 50 jam, where it is four 
times larger than for 200 pm particles. For all particle sizes, its 
increase is 50% compared to the initial value of L/d = 1. The ratio 
between convective heat flux and the heat flux into the inert par¬ 
ticle varies from 11 (cylinder, L/d = 1) to 35 (ellipsoid, L/d = 10). For 
the larger particles, this ratio becomes 0.72 (ellipsoid, L/d = 1). As 
Fig. 3 shows, qreact is influenced differently by both q inert and q con v 
for the particle shapes considered. As was mentioned before, for 
lower L/d ratios the difference between cylindrical and spherical 
(thus, elliptical particles with L/d = 1) becomes clearly visible here. 

3. Influence on experimental results 

3.1. Experimental setup and fuel 

Based on the theoretical considerations presented above, the 
model for both particle shapes has been tested on experimental 
results. A flat flame burner (a laminar flow reactor with quartz 
glass walls for optical access, see [23], Fig. 4) was used to 



Fig. 4. Drawing of the flat flame burner and the camera system. 


investigate the combustion behaviour of torrefied beech wood par¬ 
ticles. The burner, used to provide a hot oxidizing atmosphere to 
expose the fuel particles to high heating rates, was fed with meth¬ 
ane and air to achieve an atmosphere containing 7% 0 2 , 13% H 2 0 
and 7% C0 2 (N 2 for balance, volumetric concentration). The gas 
temperature along the reactor centre line was measured with a 
type-S thermo couple (corrections for radiation and gas composi¬ 
tion were applied [32]). 

In this environment torrefied beech wood particles were 
injected, which were ground and sieved to remove particles of 
more than 200 pm diameter based on the mesh size. The proxi¬ 
mate, ultimate and calorific analyses are presented in Table 3. 

Furthermore the shape of the particles was analysed before the 
combustion experiments by a Retsch Camsizer XT. This system 
uses one camera to determine the two-dimensional projection of 
the particle shape in the viewing direction, which potentially 
underestimates the L/d ratio for all particles which are inclined 
to the direction of view. Nevertheless, the fuel particles investi¬ 
gated show L/d ratios much larger than unity. The particle volume 
fraction is plotted against L/d in Fig. 5. For 30% of the sample vol¬ 
ume L/d ranges between 1 and 2 while 50% are in the L/d range 
of 2-5, and higher values complete the distribution. The data set 
used includes more than 900,000 particles for reliable statistics. 

The particle temperature and size/shape measurement was 
carried out with the imaging pyrometer system TOSCA, described 


Table 3 

Fuel analysis. 

Torrefied beech 


Prox. analysis 

Volatile matter 

(wt%, dry) 

71.03 

Fix-carbon 


28.25 

Ash content 


0.73 

Ult. analysis 

Carbon 

(wt%, dry) 

56.85 

Hydrogen 


5.51 

Nitrogen 


0.32 

Sulphur 


0.03 

Oxygen 


36.96 

Calorific value 

HHV 

(MJ/kg, dry) 

21.187 
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Fig. 5. Particle volume fraction of the initial beech wood sample. 



Fig. 6. Example of size detection in imaging pyrometry (left) and the standard 
procedure, the dark area marks the particle cross section in the line of sight 
respective it is outline (right). 


in [23]. The camera system is temperature-calibrated from 
1300-2300 K, and the spatial resolution used to measure the 
torrefied beech wood particles was 10.3 pm/px. In addition to 
the determination of particle temperature, size and shape, double 
exposed images were used to measure the particle velocity. 

Fig. 6 shows the cross section of a virtual char particle in the line 
of sight of the two-dimensional, imaging pyrometer used here. The 
diameter d sp here is calculated assuming a spherical particle shape 



Fig. 7. Gas temperature and particle velocity. 



Fig. 8. Particle temperature and diameter. 

with a projected surface equal to that of the non-spherical particle 
by counting the illuminated pixels n px and calculating 
d sp here = 2s px ^/n px /n • s px is the size calibration factor of the camera 
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Fig. 9. Scatter plots of measured spherical diameters and L/d ratios. 

































































































M Schiemann et al./Fuel 134 (2014) 53-62 


59 


system, measured with a backlight-illuminated etched grid. As is 
also shown in Fig. 6, the imaging pyrometer used here can resolve 
the particle shape directly by calculating d and L from size-cali¬ 
brated images by counting the illuminated pixels in each direction. 
This was carried out for the data presented in the following. 

3.2. Experimental results 

The gas temperature profile is depicted in Fig. 7 together with 
the mean particle velocity against the distance from the particle 
inlet. The gas temperature at the reactor centre line shows a slight 
increase directly downstream the particle inlet, followed by an 
almost linear decrease. In the relevant region (120-260 mm from 
the particle inlet), the gas temperature decreases from 1730 to 
1620 K. The mean particle velocity is given in Fig. 7 and was used 
to calculate the burning rate equations in the following. 

The particle temperatures and diameters (averages from more 
than 300 particles at each point) are presented in Fig. 8. The tem¬ 
peratures plotted against the distance to the fuel inlet show a clear 
trend from an initial temperature of 1870 K down to 1740 K. The 
standard deviation of the temperature measurement is less than 
65 K for each measurement. The particle sizes, also depicted in 
Fig. 8, are 175 pm ( d sphere ) in average, with a standard deviation 
of approx 75 pm. 

As Fig. 8 already shows, the data sets show a certain spread in 
the measured particle dimensions. Further insight into the depen¬ 
dence of the measured L/d ratio in relation to the spherical diame¬ 
ter determined is given by Fig. 9. The majority of the particles show 
L/d ratios between 1 and 2, only few particles exceed an L/d ratio 
above 3. Compared to the initial sample presented in Fig. 5, the 
fraction of high L/d ratios is reduced. This is in agreement with lit¬ 
erature describing the influence of fast pyrolysis [33], which the 
particles have experienced in the experiment presented. In general, 
smaller particles tend to larger L/d values, a fact which will be of 
importance later on. 

3.3. Burning rate analysis 

For each particle included in the data sets for Fig. 8 the burning 
rate q was calculated based on Eq. (1) and Tables 1 and 2, dividing 
the right hand side of Eq. (1) by AH. The thermal conductivity was 
calculated for the film temperature T f = (T p - T gas )/2 assuming 
free stream gas composition as described in [34]. The heat of reac¬ 
tion AH was calculated assuming that CO is the only reaction 
product. 

Fig. 10 illustrates the particle burning rates for more than 300 
particles which represent the data set collected 260 mm down¬ 
stream of the particle inlet. As the three diagrams show, the differ¬ 
ence of the three particle shapes is not visible on first sight. This 
was already predicted in Section 2.3. The still small but clearest dif¬ 
ferences can be recognized for burning rates around 10 and 1 mol/ 
(m 2 s). For q = 10 mol/(m 2 s) the burning rates assuming elliptical 
particles appear to be a few percent higher than those calculated 
for spherical particles, while those of cylindrical particles appear 
to be slightly lower. At the lower limit the burning rates of cylindri¬ 
cal particles appear to be higher than those of elliptical and spher¬ 
ical particles. However, a detailed quantitative analysis of the 
effect of particle shape assumptions for burning rate determination 
seems to be difficult based on these plots. 

To get a better impression of the effect of different L/d ratio on 
the burning rate under the assumption of different particle shapes, 
average burning rates were calculated for L/d ratios in intervals 
smaller 1.1 and from 1.1 to 2.1 in steps of 0.2, as is shown in 
Fig. 11. All particles with L/d ratios above 2.1 were summarized 
in one data point. For each L/d ratio, the averaged burning rate 
Q sphere I Qcyi respective q sphere I Q eii was calculated to point out the 


difference between the two models. Error bars indicating one stan¬ 
dard deviation of the average burning rate ratio show the scatter¬ 
ing of the data. The results calculated from all four data sets show 
roughly comparable results. The comparison between the average 
burning rates using models for spherical and cylindrical particles 
show only little deviations from unity. The burning rate calculated 
for spherical particles is only a few percent higher than that calcu¬ 
lated for cylindrical particles. However, for the data set recorded at 
the maximum distance from the fuel inlet (z = 260 mm), where the 
lowest average particle temperature is found (Fig. 8), q sp herelQcyi is 
clearly larger than 1, showing higher burning rates calculated for 
spherical particles. The ratio q sp herelQeii shows a clear tendency to 
higher burning rates determined for elliptical particles at each 
point included in the analysis, especially for larger L/d intervals. 

The results presented in Fig. 11 demonstrate the influence of 
different particle shape assumptions on the burning rates deter¬ 
mined from (2d-) pyrometric burning rate analysis. In Fig. 12, the 
L/d ratio and the burning rate ratios are grouped for d sphere - inter¬ 
vals to investigate the influence of the particle size on the burning 
rate determined for the different models. The vertical error bars 
indicate the standard deviation of each average value, while the 
horizontal bars represent one standard deviation of d sphere , which 
only acts as an indicator for the size interval presented (the whole 


spherical particles 





Fig. 10. Burning rates of particles from one measured data set under assumption of 
three different particle shapes, distance to fuel inlet 260 mm. 
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particle size distribution ( psd ) is nearly Gaussian, whereas the sep¬ 
arated intervals cannot be Gaussian as splitted parts of the whole 
psd). The average L/d ratio ranges from 1.8 to 1.4, baring a slight 
tendency to decreasing values with increasing d sphere . In contrast 
to the results presented in Fig. 11, the size-dependent burning rate 
ratios show clear tendencies. For q sp herelQcyu which was close to 
unity for size-independent L/d ratios, values of approx 1.15 are 
found for d sphere < 100 pm, while for particles with d sp here~ 
200 pm, an average value of 1-1.03 is found. The relatively large 


fraction of larger particles causes the tendency to values closer to 
unity in Fig. 11 , where only averages over L/d were used to group 
the data. For the ratio q sp herelQeiu which has already shown larger 
deviations from unity for L/d-grouped data, these tendencies are 
repeated, with q sp herelqeii < 0.9 for d sphe re< 100 pm and q spher el 
q e u < 0.84 for d sphere ~ 200 pm. For both burning rate ratios pre¬ 
sented, the standard deviations are much smaller than the differ¬ 
ence to unity, so that the effect based on the measurements is 
clearly proven. 
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Fig. 11. Ratios of burning rates determined under assumption of cylindrical and elliptical particles, each compared to the burning rate of spherical particles. 
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Fig. 12. L/d ratios and burning rate ratios for different d sphere . 
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The larger effect of the particle shape on the burning rate deter¬ 
mined for small particles points out the necessity to carefully ana¬ 
lyse the particle shape simultaneously with temperature and 
particle size. This becomes even more important, as it is well 
known, that smaller particles are preferred for the determination 
of chemical rate parameters, as small particles are less affected 
by transport limitations. Furthermore, two-dimensional ratio 
pyrometry only measures a projection of the particle shape on 
the image plane, which underestimates the real (3d) L/d ratio, 
therefore the effect would become even larger when a stereoscopic 
technique could be applied under similar conditions. 

From the results presented several questions arise which are out 
of the scope of the current work but should be discussed here 
briefly. 

First, the current paper did present char burning rates. How¬ 
ever, for the final use in CFD codes the char burning rates have 
to be converted to kinetic parameters typically modelled by an 
Arrhenius approach. The simplest assumption, frequently used in 
literature, for (oxidation only) kinetics is a single one-step Arrhe¬ 
nius expression: 

q = A 0 p n 02 /-&. (9) 

Typically, the pre-exponential factor A 0 , the activation energy E a 
and the reaction order n are calculated from burning rate measure¬ 
ments. Ways to derive the parameters for spherical particles are 
described by [9,11]. When we simplifying assume n and p 0iS as 
being constant (n typically ranges from 0.5 to 1 according to clas¬ 
sical Thiele analysis) Eq. (9) shows that an underestimation of q 
(ellipsoid particles) would lead to an underestimation of A 0 , and 
overestimation of E a . In addition it is clear from Eq. (9) that p 0iS 
is needed to calculate A 0 and E a (and n). Note that the equation 
to calculate p 0l S is itself diameter (i.e. shape) dependent [9,11]. 
Therefore, the final derivation of kinetic parameters needs further 
evaluation. 

Second, as the burning rates calculated for cylindrical and ellip¬ 
soidal particles deviate from those of spherical particles in opposite 
directions, the question is which shape should be used for a spe¬ 
cific fuel, or if one of the shapes is valid for all fuels. Coal particles 
which show even before combustion shapes similar to spheres and 
swelling during pyrolysis usually leads to further rounded particle 
shapes one could suggest an elliptical shape approximation 
because the lower limit of the elliptical shape is a sphere. For 
biomass fibres, the cylindrical shape appears to be promising. 
However, to answer the question of the “best” particle shape 
approximation, a more intense screening of different fuels with a 
wide-spread of experimental parameters (particle sizes, varying 
L/d ratios, and gas conditions) has to be carried out. Solid sampling 
at different stages of burn-out will be needed with subsequent 
microscopic analysis for further clarification. 

Third, the question arises whether the error made by assuming 
sphericity for non-spherical particles is significant compared to the 
other uncertainties implied by the evaluation procedure for the 
char burning rate and subsequent kinetic parameters. The calcula¬ 
tion of char burning rates is typically based a single-film boundary 
layer model [35 . The input parameters which are needed to solve 
the energy balance of the single film model, e.g. gas phase values 
like thermal conductivity, heat capacity, etc., particle properties 
like emissivity and char density, underlie certain uncertainties as 
well as measured values (temperature, particle shape). As the sin¬ 
gle film model, already including these uncertainties, has been 
investigated concerning its reliability to predict char consumption 
rates under pf-conditions for spherical particles [36,37] and did 
show good results compared to a detailed kinetic model, the error 
caused by particle shape negotiations (up to 18%) appears not neg¬ 
ligible. In addition, the effect of Stefan flow, although small under 


the given conditions based on the results presented in [29], has to 
be added to the model, but is not supposed to exceed the difference 
caused by particle shape effects from the current point of view. 

4. Conclusions 

The typical energy balance equation, which has been used fre¬ 
quently to determine the burning rate and thus chemical rate 
parameters for char burning kinetics, has been adapted to cylindri¬ 
cal and prolate ellipsoidal particle shapes. Theoretical calculations 
have shown that the effect of particle shape variations cannot be 
neglected when pyrometric techniques are used to characterize 
the combustion of pulverized biomass. Also, experimental data 
were analysed to determine the influence of particle shape effects 
on the burning rates of a typical torrefied wood sample. The results 
point out the following findings: 

• As the surface reaction rate of carbon q = q react /AH linearly 
depends on q react , it has to be expected that neglecting the par¬ 
ticle shape leads to a miss-interpretation of the char burning 
rate for non-spherical particles. 

• Regardless of the technique used to determine the burning rate 
of char particles, particle shape effects should be included in 
measurements used and models applied to determine kinetics 
parameters for further calculations. 

• Both theoretical and experimental results show, that the two 
particle shapes considered have opposite effects on the calcu¬ 
lated burning rate, especially for the typical particle shapes. 
The present work points out the necessity to identify which par¬ 
ticle shapes considered is most consistent with burn out mea¬ 
surements from solid samples, which then enables to use the 
optical technique to characterize the shape depending burning 
behaviour of specific fuels. 

• A very predictive technique to measure the particle shape 
simultaneously with the particle temperature could be stereo¬ 
scopic (3d) imaging pyrometry. As the results from 2d-mea- 
surements show, it is possible to measure particle shape and 
temperature simultaneously, which leads to the possibility to 
include the geometric variability and its temporal alteration in 
models to describe char burning kinetics. 

• Currently, commercial CFD-codes, which are used for predictive 
engineering of pf-boilers, do not include the particle shape in 
respect to the predicted burning rate. An inclusion of shape-cor¬ 
rected burning rate predictions for CFD-calculations should be 
developed and demonstrated. 
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